Evidence for a narrow |S| = 1 baryon state at a mass of 1528 MeV in quasi-real 

photoproduction 
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Evidence for a narrow baryon state is found in quasi-real photoproduction on a deuterium target 
through the decay channel pKg — > p7r+7r~ . A peak is observed in the pKg invariant mass spectrum 
at 1528 ± 2.6(stat) ± 2.1(syst) MeV. Depending on the background model, the naive statistical 
significance of the peak is 4-6 standard deviations and its width may be somewhat larger than the 
experimental resolution of ct = 4.3-6.2 MeV. This state may be interpreted as the predicted S=+l 
exotic 0'^{uudds) pentaquark baryon. No signal for an hypothetical B"^"*" baryon was observed in 
the pK'^ invariant mass distribution. The absence of such a signal indicates that an isotensor Q is 
excluded and an isovector O is unlikely. 
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One of the central mysteries of hadronic physics has 
been the failure to observe baryon states beyond those 
whose quantum numbers can be explained in terms of 
three quark configurations. Exotic hadrons with mani- 
festly more complex quark structures, in particular ex- 
otics consisting of five quarks, were proposed on the 
basis of quark and bag models |H in the early days of 
QCD. The hope has been that the discovery of such ob- 
jects would provide new insights into the dynamics of 
quark interactions in the strong coupling regime. Al- 
though it was hypothesized |2] that pentaquark systems 
involving heavy quarks, e.g. uudcs, offered the most 
promising prospects for isolating such exotics, experi- 
mental searches carried out at FNAL ^ found no evi- 
dence for such states. 

From quite a different point of view, it was noted 0,0 
that the Skyrme model predicts new exotic states belong- 
ing to higher SU(3) representations. Using this model, 
Praszalowicz |^ provided the first estimate of the mass 
of the lightest exotic state, M«1530MeV. Subsequently, 
an exotic baryon of spin 1/2, isospin 0, and strangeness 
was discussed as a feature of the Chiral Quark 
Soliton model ^ . In this approach 0, Q the baryons are 
rotational states of the soliton nucleon in spin and isospin 
space, and the lightest exotic baryon lies at the apex 
of an anti-decuplet with spin 1/2, which corresponds to 
the third rotational excitation in a three flavor system. 
Treating the known N(1710) resonance as a member of 
the anti-decuplet, Diakanov, Petrov, and Polyakov |^ de- 
rived a mass of 1530 MeV and a width of less than 15 MeV 
for this exotic baryon, since named the 0+. It corre- 
sponds to a uudds configuration, and decays through the 
channels — > pK^ or nK'^ . However, measurements 
of scattering on proton and deuteron targets showed 
no evidence for strange baryon resonances, and ap- 
pear to limit the width to remarkably small values of 
order an MeV. Presumably, the difficulty of experiments 
with kaon beams so close to threshold severely limited the 



sensitivity to such narrow excitations. In a recent review 
of this subject 10], an experimental search by means of 
electro-production was suggested. 

Experimental evidence for an exotic baryon first came 
recently [lH from the observation of a narrow resonance 
at 1540± lO(syst) MeV in the K~ missing mass spectrum 
for the 7n K'^K~n reaction on ^^C. The decay mode 
corresponds to a iS^-f 1 resonance and signals an exotic 
pentaquark state with quark content (uudds). Confirma- 
tion came quickly from a series of exp eriments, with the 
observation of sharp peaks [l3,ll3jb4lldlla| in the nK~^ 
and pKg invariant mass spectra near 1540 MeV, in each 
case with a width limited by the experimental resolution. 
The failure to observe a corresponding 0+''" peak in the 
pK^ invariant mass spectrum in some of these experi- 
ments 0, Q| was taken to suggest that the state is an 
isospin singlet. 

Alternative theoretical explanations have been pro- 
posed recently to explain this new exotic state. In 
one model, the Q is described as an isotensor pen- 
taquark [l7j, so that the narrow width results from the 
isospin- violating strong decay. A search for the decay of 
the isospin partners such as the Q~^~^ (uuuds) can provide 
a strong test of this idea. In a second interpretation, Kar- 
liner and Lipkin have developed a cluster model using a 
diquark-triquark configuration [l^ . in which the G"*" is 
also a positive-parity isosinglet member of an antidecu- 
plet. Thirdly, using a model based on the strong color- 
spin correlation force, Jaffe and Wilczek propose that 
the 6+ consists of two highly correlated ud pairs coupled 
to an s. In their picture the positive-parity isosinglet 
0+ lies at the apex of an nearly-ideally mixed SU(3)/ 
8f ® To7 multiplet. The narrowness of the 1530 MeV 
state is attributed to the relatively weak coupling of the 
kaon- nucleon continuum to the pentaquark [itdj^s config- 
uration. 

The baryon states at the bottom two vertices of the 
anti-decuplet must also be manifestly exotic. Strong ev- 
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idence in support of the baryon decuplet comes from the 
reported observation of an exotic 5=— 2, Q=—2 baryon 
resonance in proton-proton colHsions at -/s = 17.2 GeV 
at the CERN SPS j20|. A narrow peak at a mass of about 
1862 MeV in the S~7r~ invariant mass spectrum is pro- 
posed as a candidate for the predicted exotic S^^" baryon 

2 

with 5*=— 2, and a quark content of [dsdsu). At the 
same mass, a peak is observed that is a candidate for 
the member of this isospin quartet. The correspond- 

2 

ing anti-baryon spectra show enhancements at the same 
invariant mass. This observed mass of 1862 MeV falls 
between the predictions of Refs. 01 and 0|, although 
closer to the latter. Also, the positive parity for the O"*" 
predicted by these models contrasts with the negative 
parity suggested by the first lattice results 0, . The 
general theoretical situation is still quite unsettled. 

This Letter presents the results of a search for the 
in quasi-real photoproduction on deuterium. In addi- 
tion to corroborating some features of the state measured 
previously, the data reported here provide more restric- 
tive information related to its mass and isospin. The 
data were obtained by the HERMES experiment with 
the 27.6 GeV positron beam of the HERA storage ring 
at DESY. Stored beam currents ranged from 9 to 45 mA. 
An integrated luminosity of 250 pb~^ was collected on a 
longitudinally polarized deuterium gas target. The yields 
were summed over two spin orientations. 

The HERMES spectrometer consists of two iden- 
tical halves located above and below the positron beam 
pipe, and has an angular acceptance of ±170mrad hor- 
izontally, and ±(40 ~ 140)mrad vertically. The trigger 
was formed by either a coincidence between scintillat- 
ing hodoscopes, a preshower detector and a lead-glass 
calorimeter, or a coincidence between three scintillating 
hodoscopes and two tracking planes, requiring that at 
least one charged track appears in each of the detector 
halves of the spectrometer. 

The analysis searched for inclusive photoproduction 
of the 0+ followed by the decay 8'^ P^s ^ ptt+tt^. 
Events selected contained at least three tracks: two 
oppositely charged pions in coincidence with one pro- 
ton. Identification of charged pions and protons was 
accomplished with a Ring-Imaging Cerenkov (RICH) 
detector which provides separation of pions, kaons 
and protons over most of the kinematic acceptance of 
the spectrometer. The RICH identification efficiencies 
and cross contaminations had been determined in a 
limited kinematic domain using known particle species 
from identified resonance decays. However, because 
the RICH performance is sensitive to event topology, 
it was essential to determine these efficiencies and 
contaminations for pions and protons under conditions 
as close as possible to those of the present measurement. 
This was accomplished by means of a Monte Carlo 
simulation based on the Pythia6 generator discussed 



below. Events with the relevant topology were used to 
determine these parameters as a function of particle 
momentum. The data from the simulation indicated 
that cross contaminations would be negligible if protons 
were restricted to a momentum range of 4-9GeV/c 
and pions to a range of l-15GeV/c, the kinematic 
restrictions subsequently used in the analysis. 
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FIG. 1: Invariant mass distribution of two oppositely charged 
pions, subject to the constraints in event topology discussed 
in the text. A window corresponding to ±2 a is shown by the 
vertical lines. 

The event selection included constraints on the event 
topology to maximize the yield of the peak in the 
Mjr+TT- spectrum while minimizing its background. How- 
ever, no constraints were optimized to increase the signi- 
ficance of the signal visible in the final Mp^+^~ spectrum, 
as such optimization would have produced a spectrum 
to which standard statistical tests do not apply. Based 
on the intrinsic tracking resolution, the required event 
topology included a minimum distance of approach be- 
tween the two pion tracks less than 1 cm (the midpoint of 
which is defined as the K'g decay vertex), a minimum dis- 
tance of approach between the proton and reconstructed 
tracks less than 6 mm (the midpoint of which is de- 
fined as the production vertex), a radial distance of the 
production vertex from the positron beam axis less than 
4 mm, a z coordinate of the production vertex within the 
±20 cm long target cell of — 18cm< z < +18 cm along 
the beam direction, and a Kg decay length (separation 
of production and Kg decay vertices) greater than 7 cm. 
To suppress contamination from the A(1116) hyperon, 
events were rejected where the invariant mass Mp^- fell 
within 2 cr of the nominal A mass, where a — 2.6 MeV is 
the apparent width of the A peak observed in this exper- 
iment. 

The resulting invariant AI^+^- spectrum is shown in 
Fig. 1. The position of the Kg peak is within 1 MeV 
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of the expected value of 497.7 ± 0.03 MeV |2g. To 

search for the events were selected with a M^+^r- 
invariant mass within ±2 a about the centroid of the 
K'g peak. The resulting spectrum of the invariant mass 
of the pTT+TT" system is displayed in Fig. 2. A narrow 
peak is observed. There is no known positively charged 
strangeness-containing baryon in this mass region (other 
than the 9+) that could account for the observed peak. 
Also, the Mp^+^- spectrum corresponding to the side- 
band background adjacent to the K'g peak in Fig. 1 was 
found to be featureless. 

The non-resonant contribution to the spectrum was es- 
timated by means of a simulation using a version of the 
PYTHIA6 code tuned for HERMES kinematics 
This event generator contains no resonances in the mass 
range of Fig. 2a that decay in the channel. The 

resulting simulated spectrum is shown in Fig. 2a as the 
gray hatched histogram. The statistical precision of the 
present study is limited by the rare topology of the events 
selected. Trigger inefficiencies were not included in the 
simulation, but are believed to be small. The simulated 
spectrum falls below the data at high invariant mass 
where !]*+ resonances are known to exist |2^. There- 
fore, if Pythia6 is assumed to be capable of describ- 
ing the shape of the non-resonant contribution, it can be 
concluded that there is substantial resonant strength dis- 
tributed over the high-mass portion of the spectrum. At 
the position of the observed peak in the data, no corre- 
sponding structure appears in the simulated spectrum. 

In order to determine the centroid, width and signi- 
ficance of the peak observed in Fig. 2, three different 
models for the background were explored. For the first 
model, the PythiaB simulation is taken to represent the 
non-resonant background, and the remaining strength in 
the spectrum is attributed to a combination of known 
broad resonances and a new structure near 1.53 GeV. For 
the second model, it is assumed that the non-resonant 
background involves a large enough typical multiplicity 
that the 4-momenta of the K^g and proton are largely 
uncorrelated. In this case, this background can be sim- 
ulated by combining from different events a kaon and 
proton that satisfy the same kinematical requirements 
as the tracks taken from single events in the main anal- 
ysis. Since resonances are typically visible only as rare 
correlations between their decay particles, their contribu- 
tions will be relatively suppressed in this method. Fig. 2a 
shows that this procedure yields a shape that is very 
similar to that from the Pythia6 simulation, within the 
available statistics. By fitting a polynomial to the mixed- 
event background normalized to the PythiaB simulation, 
and then fitting this polynomial together with the am- 
plitudes of peaks for six known resonances in the 
mass range shown in Fig. 2 (dotted curves), plus all pa- 
rameters of a narrow Gaussian (dashed curve) for the 
peak of interest, a good description of the entire spec- 
trum is obtained. This procedure is intended to demon- 
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FIG. 2: Distribution in invariant mass of the p7r"'"7r~ sys- 
tem subject to various constraints described in the text. The 
experimental data are represented by the filled circles with 
statistical error bars, while the fitted smooth curves result in 
the indicated position and a width of the peak of interest. 
In panel a), the Pythia6 Monte Carlo simulation is repre- 
sented by the gray shaded histogram, the mixed-event model 
normalised to the Pythia6 simulation is represented by the 
fine-binned histogram, and the fitted curve is described in 
the text. In panel b), a fit to the data of a Gaussian plus a 
third-order polynomial is shown. 



strate that the background is consistent with known in- 
formation. The included S*+ resonances were assigned 
hxed values of M = 1480 MeV with F 55 MeV (PDG 
status = *), M = 1560MeV with F = 47MeV (**), 
M = 1580 MeV with F = 13 MeV (**), M = 1620 and 
1660 MeV with F = 100 MeV (***), and M = 1670 MeV 
with F = 60 MeV (****) Each intrinsic Brcit- 

Wigner width was taken as the midpoint of the range 
of listed measurements, and was then augmented by an 
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TABLE I: Mass values and experimental widths, with their statistical and systematic uncertainties, for the O"'" from the two fits, 
labelled by a) and b), shown in the corresponding panels of Fig. 2. Rows a') and b') are based on the same background models 
as rows a) and b) respectively, but a different mass reconstruction expression that is expected to result in better resolution. 
Also shown are the number of events in the peak Ns and the background Nt, both evaluated from the functions fitted to the 
mass distribution, and the results for the naive significance N^"/ \/ N'^" and reahstic significance Ns/SNs- The systematic 
uncertainties are common (correlated) between rows of the table. 





mass 


FWHM 






naive 


Total 


signif. 




[MeV] 


[MeV] 


in ±2(7 


in ±2(7 


signif. 


Ns ± SNs 




a) 


1527.0 ±2.3 ±2.1 


22 ± 5 ± 2 


74 


145 


6.1a 


78± 18 


4.3 a 


a') 


1527.0 ±2.5 ±2.1 


24 ± 5 ± 2 


79 


158 


6.3 a 


83 ±20 


4.2(7 


b) 


1528.0 ± 2.6 ± 2.1 


19 ± 5 ± 2 


56 


144 


4.7a 


59 ±16 


3.7(7 


b') 


1527.8 ± 3.0 ± 2.1 


20 ± 5 ± 2 


52 


155 


4.2 a 


54± 16 


3.4 a 



instrumental resolution of FWHM= 14.3 MeV added in 
quadrature. Since the E+(1580) has a width smaller than 
the instrumental resolution, it was taken to be Gaussian 
with a = 8.9 MeV. The numerical results of the fit are 
given in row a) of Table In addition, row a') shows the 
result of applying the same method to a mass spectrum 
based on an expression that is expected to provide an 
instrumental resolution improved by about 30%: 




-{Pp + Pk°)\ (1) 

where Mj^o is taken from the PDG |i2^]. The third ap- 
proach is based on the hypothesis that all of the back- 
ground strength in the observed spectrum (apart from 
the narrow peak) can be described by a slowly vary- 
ing function that extends under the feature of interest. 
Hence, the spectrum was fit with a Gaussian plus a poly- 
nomial. The appropriate degree of the polynomial used 
in the fit was determined by comparing results using 
orthonormal Chebyshev polynomials of various degrees. 
The curve shown in Fig. 2b results from a fit with a 
third-order polynomial, and rows b) and b') of Table |I| 
gives the numerical values from so fitting the two spectra 
corresponding to rows a) and a'). 

More specifically, TablelUcompares the results from the 
two fits for the centroid of the peak of interest, its width 
and the statistical significance according to two differ- 
ent prescriptions discussed below. The resulting values 
for the centroid are found to be consistent, while the 
width and significance depend on the method chosen to 
describe the remaining strength of the spectrum. Ta- 
ble ^ also lists for both fits the number of events given 
by the fitted function for the peak of interest (N^"^) as 
well as for the background {N^'^), in the invariant mass 
interval corresponding to ±2 a. The full area Ng of the 
Gaussian fitted to the peak of interest is also given with 
its uncertainty from the fit. This area itself, together 
with the width, were chosen to be explicit fit parame- 
ters to avoid the effect on the uncertainty in the area 



of correlations between the amplitude and the width or 
background parameters. All of these results are from un- 
binned maximum likelihood fits |2^ to the original event 
distributions, as it was found that the results of fitting 
the histograms shown in Fig. 2 can be sensitive to the 
choice of bin size or starting offset. 

Several alternative expressions for the significance of 
the peak observed in Fig. 2 were considered. The first 
expression is the naive estimator ^1" I \f^^ used in 
Refs. The corresponding result is 

listed in Tabled Because this statistic neglects the un- 
certainty in the background fit, it overestimates the signi- 
ficance of the peak [22| . A second estimator that was used 
in the analysis presented in Ref. HO], N"^," / yj N"^" + N^" , 
gives a somewhat lower value, but may still underesti- 
mate the background uncertainty. A third estimate of 
the significance is given by Ng/SNs, where Ng is now the 
full area of the peak from the fit and SNg is its fully corre- 
lated uncertainty. This ratio measures how far the peak 
is away from zero in units of its own standard deviation. 
All correlated uncertainties from the fit, including those 
of the background parameters, are accounted for in SNg. 
The results obtained with this expression are also given 
in Table n 

The systematic uncertainty of the mass of the state ob- 
served in Fig. 2 is estimated to be ±2.1 MeV by adding 
in quadrature a contribution of 1 MeV to account for the 
effect of using different spectrum analysis methods (cf. 
TableP) plus a contribution of 1.9 MeV from the precision 
with which the spectrometer can reproduce known par- 
ticle masses (cf. Table Hljl . The 1.9 MeV contribution to 
the systematic uncertainty accounts for both the discrep- 
ancies from the PDG |2^ mass values and the statistical 
precision of their fits. As an example, the A(1520) mass 
peak fitted with a Gaussian with free width convoluted 
with a Breit-Wigner form with its width fixed at the PDG 
value is shown in the Mpx~ spectrum of Fig. 3. The 
event selection for the spectra in this figure is the same as 
for the pK^g analysis, except the reconstructed track 
is replaced by that of the observed charged kaon, which 
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TABLE II: Masses and widths of observed invariant mass peaks for four known particles, compared to the known masses |25ll 
and the widths obtained from a Monte Carlo simulation of the spectrometer. (The widths for the H~(1321) are from mass- 
difference spectra.) The A(1520) peak of Fig. 3 was fitted with a Gaussian folded with a Breit-Wigner form whose F width was 
fixed at the PDG value (15.6 MeV). The uncertainties in the resulting mass and Gaussian a width were inflated by the factor 
\/x^ /-^dof ■ In the last row, Pcm is the momentum of each decay product in the rest frame of the decaying particle. 

K'-^-^Tv+TT- A(1116)^p7r- A{1520) ^pR- H" (1321) -» pTr^Tr" 

Observed mass [MeV] 496.8 ± 0.2 1115.70 ± 0.01 1522.7 ± 1.9 1321.5 ± 0.3 

PDG Mass [MeV] 497.67 1115.68 1519.5 ± 1.0 1321.31 ±0.13 

a Width (data) [MeV] 6.2 ± 0.2 2.6 ± 0.1 4.4 ± 3.7 3.1 ±0.3 

a Width (MC) [MeV] 5.4 2.1 3.5 2.5 

Decay Pcm [MeV/c] 206 101 244 139 (Att") 



is required to have a minimum momentum of 3 GeV. 

The mass values reported to date for the 0+ state by 
other experiments are compared to the present results in 
Fig. 4, and listed in Table ITTll . The systematic uncertain- 
ties of the DIANA and ITEP measurements were taken 
to be ±3 MeV in the absence of exphcit values quoted 
in the corresponding papers 0, 0] . By fitting the data 
with a constant, a reduced value of 12.41/6 is found, 
correspon ding to a confidence level of 0.053 as defined in 
the PDG |25j . The weighted average [2^ of the masses 
observed in all experiments is 1536.2 ± 2.6 MeV, which is 
represented by the shaded band in Fig. 4. In evaluating 
this average mass value, the quadratic sum of the statisti- 
cal and systematic uncertainties of all measurements are 
taken into account. The uncertainty of the average was 
scaled by the usual factor of square root of the reduced 



X 



Since no realistic model for the photoproduction of 




1.45 1.5 1.55 1.6 1.65 1.7 

M(pK"), M(pK^) [GeV] 

FIG. 3: Spectra of invariant mass M^j^- (top) and Mp^+ 
(bottom). A clear peak is seen for the A(1520) in the Mpj^- 
invariant mass distribution. However, no peak structure is 
seen for the hypothetical B^^ in the Mpj^+ invariant mass 
distribution near 1.53 GeV. 



exotic baryons at this experiment's energy is presently 
available, a "toy Monte Carlo" was produced to study 
the constraints imposed on the decay products by only 
kinematics and acceptance. It generates parent particles 
with specified mass and width, at vertices distributed 
according to the HERMES target gas profile. The gen- 
erated events were then passed to a full simulation of 
the spectrometer that included the performance of the 
RICH. The unknown kinematic distribution of the par- 
ent in transverse momentum Pt was taken to be Gaussian 
with a width of a = 0.4 GeV, which is typical of intrinsic 
transverse momentum of partons or transverse momen- 
tum induced by the fragmentation process, and also cor- 
responds to the root-mean-square value of the transverse 
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FIG. 4: Mass values and experimental FWHM widths ob- 
served in various experiments for the state. The inner er- 
ror bars represent the statistical uncertainties, and the outer 
error bars represent the quadratic sum of the statistical and 
systematic uncertainties. (Some uncertainties for the widths 
are not available from the other experiments.) The hatched 
area corresponds to the weighted average of all data ±1 stan- 
dard deviation. 
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TABLE III: Mass values and experimental widths for the 
state as observed in the various experiments. The present 
result is also listed. In calculating the weighted average of 
the data, the systematic uncertainties of DIANA and ITEP 
are taken to be ±3 MeV. 



T-1 .vT^ovi nrt OTi^~ 
IZjaJJcI iilitllilj 


(--1 mace 

(MeV) 


1. VV i-iiVJ. 

(MeV) 


Ref. 


SPringS 


1540 ± 10 ± 5 


25 


M 


DIANA 


1539 ±2± "few" 


9 


[12j 


CLAS (d) 


1542 ± 2 ± 5 


21 




SAPHIR 


1540 ± 4 ± 2 


25 


M 


ITEP (ly's) 


1533 ± 5 


20 


[15j 


CLAS (p) 


1555 ± 1 ± 10 


26 ±7 




HERMES 


1528 ±2.6 ±2.1 


19±5±2 




world average 


1536.2 ± 2.6 







momentum distribution of A(1116) and A hyperons that 
are observed at HERMES. The simulated acceptance was 
found to be fairly insensitive to this parameter. The lon- 
gitudinal momentum was assigned a monotonically 
falling distribution similar to what has been observed for 
A hyperons at HERMES. The acceptance for the pK^ 
final state is insensitive to this choice, but this is not the 
case for the pn^n^ final state, where drastically different 
assumptions can result in a factor of two change in the 
acceptance. The full width F of the 8^ was chosen to 
be 2 MeV, according to the limit recently derived from a 
review of KN phase shift analyses 0] . The effect of such 
an intrinsic width is small compared to that of the spec- 
trometer resolution, which was found in this simulation 
to be o- = 6.2 MeV in Mp^+^-, or 4.3 MeV in M'^o of 
Eq.[T] These values are assigned a systematic uncertainty 
of ±1 MeV, based on the level of agreement between ob- 
served and simulated widths of four known particles, as 
shown in Table ITU 

The width of the peak of interest in Fig. 2, given in 
Table is somewhat larger than the instrumental res- 
olution derived from the simulation. An attempt was 
made to repeat the fits of Table ^ using for the peak of 
interest a Breit-Wigner form convoluted with a Gaus- 
sian whose width was fixed at the simulated resolu- 
tion.^ The resulting mass values are consistent with 
those given in Table and the resulting values for the 
intrinsic width are F = 12 ± 9(stat)±3(syst) MeV in 
case a) of Table HI F = 20 ± 8(stat)±3(syst) MeV in 
case a'), F = 8 ± 8(stat)±3(syst) MeV in case b), and 
F = 13±9(stat)±3(syst)MeVincaseb')- The systematic 
uncertainties here correspond only to the 5a = ±lMeV 
uncertainty in the instrumental resolution, which was dis- 



^ The software tool RooVoigtian was used, for unbinned fitting. 



cussed above. 

In view of the speculation that the observed resonance 
is isotensor 0, the possibility that the 6++ partner is 
present in the Mpx+ spectrum was explored. Although 
Fig. 3 shows a clear peak for the A(1520) in the Mpx~ 
invariant mass spectrum, there is no peak structure ob- 
served in the Mpx+ invariant mass distribution. From 
a fit (curve in Fig. 3) of the Mpx+ distribution using a 
free polynomial plus a Gaussian with the fixed location 
(±5 MeV) from the fit in Fig. 2b and a fixed width corre- 
sponding to the simulated peak width of cr = 4.5 MeV in 
this decay channel, the Gaussian area for a hypothetical 
0++ peak is found to be —40 ± 30 events. This result is 
robust against varying the order of the polynomial back- 
ground. It corresponds to an upper limit of zero counts 
at the 91% confidence level. 

The failure to observe a 6++ suggests that the 9 is 
likely to be isoscalar. However, in the situation more 
probable at lower beam energy that the Q is produced 
only via the exclusive reaction 7 ± p ^ Q + K without 
any other hadrons in the final state, the following limi- 
tations would apply to deductions about its isospin. Un- 
der the assumption of isospin symmetry, selection rules 
limit the transition amplitude for forming a tensor Q to 
a single reduced matrix element for an isovector transi- 
tion. In this case, production of the 8++ and the 6+ 
are expected to have comparable strength, and the fail- 
ure to observe the 8++ rules out the 1 = 2 assignment. 
Production of an isovector 8 would arise from a sum of 
three reduced amplitudes with unknown magnitudes and 
phases. With only model-dependent values for these am- 
plitudes, no precise statement can be made about the 
relative yields of the 8++ and 8+. However, because a 
nearly complete cancellation is improbable, the failure to 
observe the 8++ indicates that an isovector 8 is unlikely. 

Estimates of the spectrometer acceptance times effi- 
ciency from the toy Monte Carlo simulation can be used 
to estimate some cross sections. Using the assumptions 
about the initial kinematic distribution described above 
and assuming that the decay angle distribution is flat, 
these acceptances were estimated to be 1.5% for both 
A(1520) -> pK~ and 8++ ^ pK+, and 0.05% for 
8+ ^ pf^s- Taking the branching fraction of the 8+ 
to pKg to be (1/2) • (1/2) (to account for competition 
with both the nK^ channel and K^), the cross section 
for its photoproduction is found to range from about 100 
to 220 nb ±25%(stat), depending on the model for the 
background and the functional form fitted to the peak. 
The cross section for photoproduction of the A(1520) is 
found to be 62 ± ll(stat) nb. Hence the ratio of the 8+ 
cross section to that for the A(1520) is found to lie be- 
tween 1.6 and 3.5. All of these estimates are subject to 
an additional factor of two uncertainty, to account for 
the assumptions about the kinematic distribution of the 
parents used in the simulation as explained above, and 
neglected trigger inefficiencies. 
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In conclusion, evidence has been obtained in quasi- 
real photoproduction on a deuterium target for a nar- 
row baryon state in the pKg invariant mass spectrum 
at 1528 ± 2.6(stat) ± 2.1(syst) MeV. Depending on the 
background model, the width of the observed peak may 
be larger than the experimental resolution of a =4.3- 
6.2 MeV. Fitting the peak with a convolution of a Breit- 
Wigner shape with a Gaussian representing the simu- 
lated instrumental resolution yields an extracted intrinsic 
width r = 17 ± 9(stat)±3(syst) MeV (the average of the 
results from cases a' and b' of Table QJ. The significance 
of the observed state expressed as N^'^/^/N^ ranges 
from 4.2 a to 6.3 cr, and expressed as Ns/SNg ranges from 
3.4 a to 4.3 a, again depending on the model for the back- 
ground. This observation provides further evidence for 
the existence of a narrow baryon state with jS"! = 1 and 
a mass in the region where such a feature was observed by 
earlier experiments. Formally, the difference between the 
value for the mass derived here and that from the other 
experiments reduces the confidence level of the combined 
fit to all mass data from 0.57 to 0.053, a value still typical 
for well-established particles [2^. There is no identified 
S*+ state with 5*=— 1 in the invariant mass region be- 
tween 1500 and 1550 MeV. Therefore, the state observed 
here may be interpreted as the predicted exotic 0"*" pen- 
taquark S'=-|-l baryon. The absence of a corresponding 
signal in the pK~^ invariant mass spectrum indicates that 
the observed G'^ is not isotensor and is probably an isos- 
inglet. 
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